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Abstract 

We present a phenomenological approach to dispersion in nonlinear 
elasticity. A simple, thermomechanically sound, constitutive model is 
proposed to describe the (non-dissipative) properties of a hyperelastic 
dispersive solid, without recourse to a microstructure or a special ge- 
ometry. As a result, nonlinear and dispersive waves can travel in the 
bulk of such solids, and special waves emerge, some classic (periodic 
waves or pulse solitary waves of infinite extend), some exotic (kink or 
pulse waves of compact support). We show that for incompressible 
dispersive power-law solids and forth-order elasticity solids, solitary 
waves can however only exist in the case of linear transverse polar- 
ization. We also study the influence of pre-stretch and hardening. 
We provide links with other (quasi-continuum, asymptotic) theories; 
in particular, an appropriate asymptotic multiscale expansion special- 
izes our exact equations of motion to the vectorial MKdV equation, 
for any hyperelastic material. 
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1 Introduction 

The interplay of short and long range interactions in a physical phenomenon 
leads to the existence of coherent structures that may play an important 
role in the determination of the dynamical and thermodynamical properties 
of real materials. That is why highly sophisticated nonlinear lattice models 
taking account of long-range interactions have been extensively used to in- 
vestigate the complex behaviour of many materials, from crystalline solids 
to rubber-like elastomers. These models lead to major achievements in the 
study of those phenomena where inherent characteristic lengths are a funda- 
mental ingredient of the mechanical behaviour, such as fracture mechanics 
or nano-structures (carbon nanotubes or bio-polymers). 

The mathematical models of atomic lattices involve large nonlinear sys- 
tems of ordinary differential equations and their treatment requires an enor- 
mous computational effort. Moreover, it is hard to conceive that such equa- 
tions may be amenable to a simple, yet detailed, mathematical analysis. On 
the other hand, standard continuum theories - which assume that the lattice 
parameter is equal to zero - are more appealing than their discrete counter- 
parts because they provide a powerful synthetic tool of analysis, description, 
and prediction. Their weakness is that they become all but useless at short 
wavelengths (a well-known example of this flaw is recalled below.) Hence 
the quest for quasi- continuum models, which generalize standard continuum 
mechanics to incorporate characteristic material lengths. 

A well-known example of quasi-continuum model is the so-called strain- 
gradient elasticity theory, which incorporates higher derivatives of the dis- 
placement. This theory is clearly connected with the higher terms in the 
naive (Taylor) expansion in terms of small h of a lattice. To reveal this con- 
nection, let us consider the vibration of a single 1-D lattice, via the following 
Hamiltonian 



where y n is the displacement from equilibrium of the n— th particle, h is the 
equilibrium inter-particle distance, and P is the inter-particle potential. The 
equations of motion derived from (II. lft are 



where u n = (y n — y n -i)/h and T(-) = P'(-)/h. The continuum limit of this 




n 






equation, up to 0(h 4 ), is 



u tt = [T(u)] 



XX 





2 



Here it is the last term (fourth-derivative with respect to space) which is 
usually hailed as being the first correction needed to account for the dispersive 
effects due to discreteness. 

Mindlin [1] introduced phenomenological strain-gradient theories of con- 
tinuum mechanics in the early 1960s; Green and Rivlin [2] extended them to 
include strain-gradients of any order; and Toupin and Gazis [5] established a 
rigorous correspondence between strain-gradient theories and atomic lattices. 
Equation (11.31) presents an example of this connection. This equation is actu- 
ally quite popular in all those applications of continuum mechanics where in- 
herent characteristic lengths of materials may not be neglected, even though 
it is associated with several unpleasant problems. For instance, because of 
the presence of a fourth-order space derivative term, additional boundary 
conditions are required. Also, the initial-value problem for (II .3p is ill-posed. 
Indeed, the associated dispersion relation is 

co 2 „o h? 

12" 

in the linear case, where c is the speed of sound, u is the frequency, and K 
is the wave number. Clearly here, instability develops at short wavelengths 
(when K 2 h 2 > 12). Rosenau [I] overcame this problem by starting from the 
exact dispersion relation for the linear discrete system (11.21) . that is 

uj 2 4 

- = -sin 2 (/^/2), (1.5) 

to deduce a regularizing expansion in the parameter h, 

uj 2 K 2 

;i.6) 



- = K 2 - (1.4) 



c 2 1 + h 2 K 2 /12 

It is clear that equations (II. 4p and (II. 6p are equivalent at large wavelengths 
(small K), but not at short wavelengths (large K), where the latter remains 
bounded. Now rewrite the continuum limit of (II. 2p as 

u tt = L A D 2 [T(u)}, where D = d x , L A = l + ^D 2 + ^-D 4 + . . . (1.7) 

12 3o0 

Then observe that L A is an invertible Schroedinger operator, with 

h 2 h 4 

L~} = 1 — —D 2 + ——D A + (1.8) 

A 12 240 K J 

Applying the operator L^ 1 to both sides of (II. 7p gives, up to 0(h 4 ), the 
equation 

h 2 

u tt = [T{u)} xx + — u xxtt . (1.9) 
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In the linear limit, this equation possesses the (bounded) dispersion relation 
(11. 6p . Also, the Cauchy problem for this equation is well-posed and we do 
not require additional boundary conditions [9]. 

Rubin, Rosenau, and Gottlieb [10] elaborated a phenomenological 3D 
counterpart to the regularization procedure of Rosenau, by modifying the 
Cauchy stress tensor T and the free energy if) to include dispersive effects, 
while keeping the other thermomechanical entities (the entropy, the internal 
production of entropy, and the entropy flux) unchanged. Destrade and Sac- 
comandi [7], [8] recently showed that the material model of Rubin, Rosenau, 
and Gottlieb falls within the theory of simple materials, which is a theory 
where the Cauchy stress tensor depends only on the history of the deforma- 
tion gradient. The consequence is that it is possible to introduce inherent 
characteristic lengths without having to introduce exotic concepts such as 
the concept of hyperstress, necessary to deal with strain-gradient theories. 

The aim of the present paper is to provide a general treatment for the 
propagation of large amplitude transverse bulk waves in the framework of a 
phenomenological theory. Hence we aim at generalizing the celebrated study 
of Gorbacheva and Ostrovsky |5j, which presented results for the continuum 
limit of a one-dimensional lattice with elastic bonds under longitudinal stress. 
Their results were obtained only for a special law of the elastic bond and for 
a single progressive wave (in a continuum KdV limit of the lattice). Here we 
consider a general constitutive law for the strain-energy of the material and 
we also consider the possibility of nonlinear dispersion. Thus we expect to be 
able to extend the results of Gorbacheva and Ostrovsky to Murnaghan's ma- 
terials, to Ogden's materials, and to many other popular phenomenological 
response functions. Focusing on incompressible solids, and using the meth- 
ods developed by Destrade and Saccomandi [B] [8], we present the governing 
equations in Section 2 and derive in Section 3 a single complex wave equation 
for finite-amplitude transverse principal waves, valid for any dispersive solid. 
In Section 4 we specialize the equation to the case of a strain-hardening 
power-law solid, and discuss the influence of polarization, dispersion, pre- 
stretch, and strain-hardening. We find some exotic solutions such as pulses 
or kinks with infinite or even compact support. However the analysis soon 
reveals that these localized solutions are the exception rather than the rule. 
The analysis is exact and does not rely on regularized expansions. When we 
do perform such an expansion (Section 5) we recover, within an appropriate 
asymptotic limit, a "vector Modified Korteveg-deVries equation" (MKdV). 
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2 Governing equations 



Let the motion of a body be described by x = x(X, t), where x denotes the 
current coordinates of a point occupied at time t by the particle which was 
at X in the reference configuration. The associated deformation gradient is 
F(X, t) = <9x/<9X and the spatial velocity gradient is L(X, t) = <9v/<9x, where 
v = dx/dt is the velocity vector. We consider a material whose mechanical 
behaviour is described by a given Cauchy stress tensor, T E say. Standard 
continuum theories (such as the linear theory of elasticity or the Navier- 
Stokes equations of fluid mechanics) are intrinsically size-independent; to 
overcome this shortcoming we modify the standard Cauchy stress tensor T E 
so that the equation for the full Cauchy stress tensor is 

T = T E + T D , (2.1) 

where the new stress tensor term, T D , is introduced to take into account 
dispersive effects. Guided by preliminary work [6], [8], we take it in the form 

T D = «(D-D)[A 2 - Aj], (2.2) 

where D is the stretching tensor, A x and A 2 are the first two Rivlin-Ericksen 
tensors, 

D = (L + L T ) /2, Ai = 2D, A 2 = A x + A : L + L T Ai, (2.3) 

and the dispersion material function a = a(D • D), must be positive due to 
thermodynamics restrictions. 

It turns out that the term (12. 2p is exactly the one proposed by Rubin et 
al. in [TP] , and that T D is a special case of the extra Cauchy stress tensor 
associated with a non-Newtonian fluid of second grade, which is 

uA l + a 1 A 2 + a 2 A 2 l , (2.4) 

in general, where v is the classical viscosity and a±, a 2 are the microstructural 
coefficients. Note that f)2.2p is of the same form as (12 .4p when v = and cti + 
a 2 = 0; here the first equality means that the material is non-dissipative, and 
the second makes the model compatible with the laws of thermodynamics, see 
Fosdick and Yu [Hj for details (If we let v ^ 0, then we have the possibility 
to include dissipation as in the classical Navier-Stokes theory, but we do not 
pursue that alley here.) The coincidence between (I2.2p and the model of 
Rubin et al. is completed once we identify the dispersion material function 
a with the derivative of the Helmholz free energy ip which they introduced 
to model dispersion: 

a = V'(D-D). (2.5) 



5 



Now we restrict our attention to homogeneous, isotropic, compressible 
elastic solids. The response of such materials from an undeformed reference 
configuration is described by the constitutive relation 

T E = + v^— I + -7=7^ + 2 v^— B-\ 2.6 

\VhdI 2 dI 3 J y/hdh dh 

where I is the identity tensor, B is the left Cauchy-Green strain tensor defined 
by B = FF T , and the strain-energy function E is a function of the first three 
invariants I\, 1%, I3 of B, defined in turn by 

J 1= trB, J 3 = detB, J 2 = \ [If - tr (B 2 )] = I 3 tr (B" 1 ). (2.7) 

When the material is incompressible, det F = 1 at all times (every motion 
is isochoric), and we replace ( 12. 6p with 

T E = -pi + 2SiB + 2E 2 B -1 , (2.8) 

where now E = E(Ii, I 2 ) only (because I3 — 1 at all times), p is the Lagrange 
multiplier due to the incompressibility constraint, and E, = dH/dli. 

For the rest of the paper, we focus on incompressible solids and we try 
to work in all generality for the elastic part of the Cauchy stress tensor. 
For applications and illustrative examples, we shall consider for instance the 
power-law model proposed by Knowles 



2b 



1 + - (h - 3) 
n 



(2.9) 



where \l > is the infinitesimal shear modulus, and b > 0, n > are material 
constants. The neo-Hookean model, for which the generalized shear modulus 
is constant, corresponds to the case n = 1. The material modelled by (12 .9p 
is softening in simple shear when n < 1 and it is hardening in simple shear 
when n > 1; hence Raghavan and Vorp [16] recently modelled the elasticity 
of aortic abdominal aneurysms by taking n = 2; also, on letting n — > 00 
we find Fung's exponential strain-energy density, widely used to describe 
soft biological tissues. We also mention models which are popular in the 
so-called weakly nonlinear elasticity theory, based on truncated polynomial 
expansions of the strain energy density. For instance, finite element packages 
often propose a polynomial strain energy potential; the second degree version 
of this potential is 

E = C 1 o(/i-3) + Coi(/2-3) + C 20 (/i-3) 2 + Co2(/2-3) 2 + C 11 (/ 1 -3)(/ 2 -3), 

(2.10) 
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where the Cy are five material parameters to be adjusted for optimal curve 
fitting. Note that this strain energy density coincides with the power-law 
n = 2 model when I\ = I 2 . Another important model comes from the 
nonlinear acoustic literature, where the strain energy function is expanded 
up to the fourth order in the strain in order to reveal nonlinear shear waves. 
In that framework, Murnaghan's expansion [T7| is often used: 

A + 2/i , I + 2m , 42 2 

£ = — - — ii - 2//z 2 H — i 1 - 2mi\% 2 + m 3 + v\i x + v 2 i x i 2 + v&ih + v±i 2 , 

(2.11) 

where A, p are the Lame moduli, I, m, n are the third-order moduli, and Ui, 
v 2 , V3, ^4 are the fourth-order moduli. Here we used the first three principal 
invariants i±, i 2 , 13 of E, the Green-Lagrange strain tensor; they are related 
to the first three principal invariants I\, I 2 , I3 by the relations, 

J 1 = 2z 1 + 3, I 2 = 4ii + 4i 2 + 3, h = 2ii + 4i 2 + Si 3 + 1. (2.12) 

When the material is incompressible, i.e. the relative portion of energy stored 
in compression is negligible, Hamilton et al. [18] proposed the following 
reduced version of the expansion (12. lip . 

£ = —2fii 2 + ni 3 + vj^. (2-13) 

Similarly, in the examples we specialize the dispersion material function 
of (I2.2p and (12. 5 p to the following simple form, 

a = pPo + pP 1 (D-T>), so that ^ = p/3 (D ■ D) + P y (D ■ D) 2 , (2.14) 

where /3q and @i are positive material constants. 

Remark 1: We emphasize that although the stretching tensor D and 
its objective time derivative appear explicitly in the expression (I2.2p . the 
dispersive stress tensor T D nonetheless does not contribute to dissipation. 
In that respect, its action is analogue to a gyroscopic force term in classical 
mechanics. 

Remark 2: That T D gives a constitutive equation compatible with the 
class of simple materials is a happenstance. If we were to consider further 
terms in the approximation of the discrete linear dispersion relation (11.51) . 
this property would no longer be true. Hence, to approximate (11.51) to order 
0(h 6 ), it is necessary to introduce the gradient of the strain. Then, the quasi- 
continuum modelling of the material has to turn to more complex theories 
and the mathematical simplicity and feasibility of the present model is lost. 
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Remark 3: We note that the regularized equation ( II. 9p of Rosenau, and 
the constitutive model (12.1 p both involve derivatives with respect to time, 
and are thus mainly suitable to study weak-nonlocality in the framework of 
wave propagation. 

Remark 4: We focus mostly on incompressible solids because it is a 
convenient way to bypass the problem of a coupling between transverse and 
longitudinal waves, and to build a general three-dimensional generalization 
of the one-dimensional lattice theory used by Gorbatcheva and Ostrosky [5], 
where transverse waves were considered coupled to a longitudinal static de- 
formation (note that Cadet [11], [13] considered a lattice model with longitu- 
dinal/transverse wave coupling.) In compressible materials, transverse waves 
are in general coupled to the longitudinal wave, except for those special forms 
of the strain energy density function which satisfy X1 + 2X2 + S3 = const., see 
the survey [19J for details. The fourth-order elasticity model (12.131) clearly 
fits into that category. 



3 Finite amplitude waves in incompressible 
solids 

In this Section we consider the following class of motions 

1 1 
x = \~2X + u(z,t), y = \~2Y + v(z,t), z = XZ, (3.1) 

where A is the pre-stretch in the Z direction. These are motions describing 

a transverse wave, polarized in the {XY) plane, and propagating in the Z 

direction of a solid subject to a pure homogeneous equi-biaxial pre-stretch 

along the X, Y, and Z axes, with corresponding constant principal stretch 
_i 1 

ratios A 2 ; A 2, and A, respectively. Here u and v are yet unknown scalar 
functions of z and t. 



3.1 Equation of motion 

The geometrical quantities of interest 
B and its inverse B _1 , given by 

A -1 + X 2 u 2 z X 2 u z v z X 2 u z 
X 2 u z v z A~ x + X 2 v 2 X 2 v z , 
X 2 u z X 2 v z X 2 



are the left Cauchy- Green strain tensor 



A — Xu z 

A -Xv z 
-Xu z -Xv z X(u 2 z + v 2 ) + X- 2 

(3.2) 
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respectively, where the subscript denotes partial differentiation. The first 
two invariants of strain are thus 



h = 2A- 1 + A 2 + \\u\ + v 2 z ), I 2 = 2X + X- 2 + \{u 2 z + v 2 z ), (3.3) 

and they are clearly functions of u 2 + v 2 only, which itself depends on z and 
t only. The kinematical quantities of interest are A 1; A 2 , and A 2 , given by 
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Vztt 


2« + v 2 zt ) 



(3-4) 

respectively. Hence we find that D-D, and thus a and ip, are functions of 

u lt + v zt on iy : 

n-n = l(u 2 zt + v 2 zt ). (3.5) 

Now the equations of motion, in the absence of body forces, are given 
in current form as: div T = pd 2 x./dt 2 , where p is the mass density. For an 
incompressible material, they read 



dp | dT 13 
dx dz 



putt, 



dp | dT 23 
dy dz 



PVtt, 



dT- 



33 



dz 



0. 



(3.6) 



Differentiating these equations with respect to x, we find p xx = p yx = 
Pzx = 0, so that p x = q\(t), say. Similarly, by differentiating the equations 
with respect to y, we find p y = (feW, say. Now the first two equations reduce 
to 



-qi{t) + {Qu z ) z + (au zt t) z = pu tt , 



-q2(t) + (Qv z ) z +(av ztt ) = pv tt , (3.7) 



and the third equation determines p. Here, Q = Q(u 2 + v 2 ) is the generalized 
shear modulus, defined by 



Q = 2(A 2 £ 1 + AE 2 ) 



(31 



and a = a{u 2 zt + v 2 t ) is the material function describing dispersive effects. 
From the generalized shear modulus, it is possible to derive the shear stress- 
shear strain law for the motion of an incompressible material. Strong exper- 
imental evidence — at least for rubber-like materials — suggests that Q > 
for any isochoric deformation [20J. 

Following Destrade and Saccomandi [6] , we take the derivative of Eqs. (13.71) 
with respect z, we introduce the strains U = u z , V = v z and the complex 
function W = U + iV, and we recast Eq. (13. 7p as the following single complex 
equation, 

(QW) ZZ + {aW tt ) zz = pW tt , (3.9) 
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where Q is now a function of U 2 + V 2 alone: Q = Q{U 2 + V 2 ), and a is now 
a function of U 2 + V 2 alone: a = a(U 2 + V 2 ). 

It is convenient to decompose the complex function W into its modulus 
VL (say) and its argument 9 (say), 

W(z,t) = n(z,t)e w{z ' t \ so that Q = Q(tt 2 ), a = a(n 2 + tt 2 9 2 ). (3.10) 

Here Q is the strain wave amplitude and 9 is the ang/e o/ polarization. Note 
that the combinations of (13.81) and (13.31) and of (12. 5p and (13. 5p give the 
relations: 

dS d?/' 

g = 2 d(^' a = 2 dWTW (3 ' n) 

3.2 Travelling transverse waves 

Now we focus on travelling wave solutions in the form 

W = W(s) = n(s)e id{s \ where s = z - ct, (3.12) 

c being the speed. This ansatz reduces (\3M) to (QW)" + c 2 (aW")" = pc 2 W". 
We integrate it twice, taking each integration constant to be zero in order to 
eliminate the rigid and the homogeneous motions. We end up with a vector 
nonlinear oscillator, 

(Q - pc 2 ) W + c 2 aW" = 0. (3.13) 
Then, separating the real part of this equation from the imaginary part gives 

[Q(n 2 ) - pc 2 ] n + a (c 2 ti 2 + C W 2 ) c 2 (n" - 9' 2 n) = o, 

(n 2 9')' = 0. (3.14) 

We integrate the latter equation to 9' = IQ~ 2 , where J is a constant, and we 
substitute into the former equation to get 

[Q(Q 2 ) - pc 2 ] Q + a (c 2 ti 2 + c 2 / 2 ^- 2 ) c 2 (Q" - I 2 Q-' 3 ) = 0. (3.15) 

Multiplying across by fl', and using the connections (13.111) . we find the energy 
first integral, 

S(fi 2 ) - \pc 2 n 2 + V (c 2 n' 2 + c 2 I 2 n- 2 ) = E, (3.16) 
where E is a constant. 
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Now a standard analysis using phase-plane theory [5] or a Weierstrass 
"potential" theory [21], [22] gives a classification of all possible travelling 
waves solutions. Note that the choice 1 = gives 9 = const., and then we 
obtain the special case of a linearly-polarized transverse wave. 

Eq.( l3.16|) is the central equation of this paper. It is valid for any dispersive 
incompressible solid with constitutive equations (12. ip . (I2.2p . ( I2.5p . and (12.81) . 
It is exact and no approximations nor regular expansions have been made to 
derive it. Before we move on to a constitutive example, we check that the 
equation is coherent with quasi- continuum theories. 



3.3 Link with nonlinear lattice theory 

The constitutive model (I2.ip . (12.21) . (12. 8p is grounded on quasi-continuum 
models, and we expect that (13.161) covers atomistic theories. In fact the 
correlation is remarkably simple: it suffices to take a = const, (i.e. to take 
fix = in (I2.14p ) and to take Q = /i + pifl 2 where /i , p\ are positive 
constants (i.e. £ is of the polynomial form ( I2.10p . see [S]). Then Eq. (l3.15p 
reduces to 



-a 1 - ^ i - 



r'2 



tt = 0. (3.17) 



_pc 2 /3 fa \ pc 2 J n\ 

Gorbatcheva and Ostrosky [5j and Cadet [H] obtained this equation for a 
mono-atomic chain with interatomic interactions by a power series expansion 
up to the fourth order over the lattice parameter. 

4 Travelling waves in hardening power-law 
solids 

To investigate the effect of nonlinearity in the elasticity and in the disper- 
sion, we first take a look at the situation for transverse travelling waves in 
dispersive hardening power-law solids, for which E is given by ( 12. 9 p and i/j 
by (I2.14p . We start with the material at n = 2, and we indicate the trends 
for n > 2. 

4.1 Linearly polarized transverse waves 

When n = 2 and the material is not pre-stretched (A = 1), we find that 
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Then, the motion of linearly-polarized waves (J = 0) is governed by the 
following specialization of ( I3.16p . 

(ji - P c 2 + ^-n 2 )n 2 + P c 2 (f3 n' 2 + ^n' 4 ) = 2E. (4.2) 

First we investigate the situation when /3q ^ 0, fix — 0. Then we recover 
the classic solutions of quasi-continuum modelling of nonlinear lattices. To 
show this, we call Qq the value of Q when Q 2 = 0, and we favour Q i n stead 
of E as an arbitrary constant (they are related through 2E = (/i — pc 2 )f2p + 
fibQ^/A). Then we perform the following change of variable and change of 
function, 

and the non-dimensional version of (13.161) at /?i = is 



u* = -(u 2 - 1) 



(4.4) 



For the choice E — (or equivalently fi' = at = 0), the changes of 
variable and of function are 




and the non-dimensional equation is 

J 2 = -uj 2 (lo 2 - 1). (4.6) 

In the latter case, the wave speed is necessarily supersonic with respect to 
an infinitesimal bulk wave (i.e. pc 2 > p) but otherwise arbitrary, and the 
solution is a well-known pulse solitary wave, o>(£) = sech £, shown on Figure 
1(a). In the former case (I4.4p . the solutions are periodic waves in general. 
When the arbitrary speed is such that pc 2 < p+pbQl/4, the amplitude varies 
between -1 and +1 and the solution is found explicitely in terms of a Jacobi 
elliptic function as 

w(0 = cn (Vl-7 2 £ I l/Vl-7 2 ) , where 7 = ^ 1 + 4(p - pc 2 ) / (pbQ 2 ) . 

(4.7) 

Notice that here the finite amplitude periodic wave can travel at the same 
speed as that of an infinitesimal shear bulk wave yfpTfp. Figure 1(b) displays 
this wave in solid curve when 7 = 1 (and then pc 2 = p and the wavelength 
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is ~ 5.24) and in dashed curve when 7 = 2 (and then the wavelength is 
~ 2.97). When the speed is such that pc 2 > p+pbQl/4, there are two possible 
waves, one with amplitude varying between — a/ 4(pc 2 — p) / (pbVll) ~~ 1 an -d 
-1, the other with amplitude varying between a/ 4(pc 2 — p)/ (pbQl) — 1 and 
+ 1. Finally, there is a special wave for the choice pc 2 = p + p&f2Q/4, because 
then Eq. fl4.4p has the pulse solitary wave solution: oj = sech £. 

A different picture emerges when /3 = 0, Pi 7^ 0. In the case E 7^ 0, the 
changes of variable and of function, 



pb 



1/4 



O(s) 



(4- 



give the non-dimensional governing equation 



00 



'4 



(a; 2 - 1} 



w 2 + l 



fibfl^ 



(pc 2 - p) 



whereas the choice E = and the changes 



pb 

pWi 



1/4 



«(0 



p6 



[4(pc 2 -p)J 



1/2 



n(s) 



(4.9) 



(4.10) 



give 



2/ 2 

-a; (u; 



1). (4.11) 

For (14. 9p . the solutions are periodic waves in general, just as when (3q ^ 0, 
0i = 0, see above. Equation (14.111) however has supersonic solitary wave 
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-1 a 1 2 3 -1 u 1 2 3 4 5 

(a) kink (b) pulse 



Figure 2: Compact solitary waves in a power-law dispersive solid. For the 
kink, the solution cu(£) is defined on the compact support [0, 7T/4 1 / 4 ]; for the 
pulse, on [0,2tc/4V% 

solutions with compact support, a rare occurrence for waves in solids. Indeed, 
fl4.1ip can be integrated formally to give [8], = X _1 (^), where X~ l is 
the inverse of a monotonic function defined in terms of a hypergeometric 
function as X{x) = 2(a; 2 ) 1//4 Hyp 2 Fi[l/4, 1/4, 5/4, x 2 ]. Figure 2 shows how we 
can construct (weak) solutions with finite support measures as a compact 
solitary kink wave or as a compact solitary pulse wave. Notice also that these 
solutions also occurs for (14. 9p at the special speed given by pc 2 = + fibQl/4. 

When fa ^ 0, j5\ ^ 0, the governing equation is a quadratic in (fl') 2 . Its 
resolution is straightforward and the results are similar |22j to those in the 
case where (3q ^ 0, f3\ = that is, solitary and periodic waves on infinite sup- 
port measures (no compact waves). In this case we checked, using the same 
arguments as those used in [23], that when the ratio fio/ fi\ approaches zero 
the tails of the non-compact localized wave decay more and more rapidly and 
the solution approaches more and more the limiting compact wave (obtained 
at p = 0). 

4.2 Plane polarized transverse waves 

Now what happens when the wave is not linearly polarized (J ^ 0)? Then 
we find that the presence of the I 2 Q~ 2 term in the equations upsets the 
delicate balance between nonlinearity and dispersion which allowed for the 
appearance of localized, and even compact, solitary waves. Only periodic 
solutions exist when 1^0. To show this, we take the case (3 ^ 0, j3x — 0, 
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A = 1. 

When E ^ 0, the changes of variable and function (14.31) give the following 
non-dimensional equation, 



UJ 



uj- 



uj a + 



1-4 



pc — /i 
pbVt^ 



u 2 -A 



pc 



.2 r2 



pbVLl 



(4.12) 



Now recall that the existence of a localized (pulse or kink) solitary wave is 
subordinated to the existence of a double root in the right hand-side of these 
equations [21], [22J. According to the choice for c, Eq.f l4.12p can be written 
in the form 



UJ 



UJ 



(w 4 ± 7 \; 2 -5 4 ), 7 = 



1-4 



pc — p 



pbVt 



1/2 



5 = 



pc 



.2 t2 



1/4 



pbtt b 

'•13) 

where the plus sign is taken when pc 2 < 1 — 4(pc 2 — p)j (phQ^) and the minus 
sign is taken when pc 2 > 4(pc 2 — p)/(pbQl) — 1. A simple analysis of these 
equations reveals that their right hand-side cannot have a double root and 
remain positive at the same time. It follows that here, there are no solitary 
waves, only periodic waves. 

Similarly when E = 0, the changes of variable and function (14. 5p give an 
equation which can be rewritten as 

2 



UJ 



'2 



(w 6 - uj 4 + 5) , where 5 = pc 2 (3 I 2 



pb 



, • • (4-14) 

W ' [A(pc 2 -p)\ 

Here the discriminant of the cubic in uj 2 is —5(275 — 4); it is zero when 
275 = 4, but the corresponding right hand-side of f)4.14p is negative. 

We conducted the same analysis when /?i ^ and checked that again, 
the other solitary localized waves (pulses, kinks, compact-like) disappear due 
to the introduction of the singular term I 2 /Q 2 . A direct and straightforward 
analysis of ( I3.16P makes it clear that this conclusion does not depend on the 
particular choice of £ made here. When the generalized shear modulus sat- 
isfies the empirical inequality Q > 0, localized waves are bound to disappear 
when 1^0. 



4.3 Pre-stretch 



When the solid is pre-stretched, the principal strain invariant I\ is given by 
(13. 3p and equation (I3.16P is changed accordingly. In the case 1 = 0, Eq.f l4.2p 
is replaced with 



(3 



pc 2 (t3 n' 2 + '^n'*) = p\ 



l+ U -(2X^ + X 2 ) 



(n 2 



n 2 Q ) + p-\ 4 (tf-n 4 0J , 



(4.15) 
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where f2 is the value of f2 when VL 2 = (it is related to E through 2E = 

mi)-) 

Hence the analysis is hardly modified by the introduction of pre-stretch. 
Methodologically and qualitatively, the results of Sections 4.1 and 4.2 ap- 
ply here; quantitatively, it makes no sense to compare the unpre-stretched 
situation with the pre-stretched case because fl and c are arbitrary. 



4.4 Other strain-hardening power-law solids and fourth- 
order elasticity 

When n in (12. 9p is an integer other than 2, the analysis is not overly modified. 
For instance, the factorizations occurring at n = 2 on the right hand-side of 
(14. 4p and (14. 6 p are still in force, and the bracketed term is now a polynomial 
of degree n — 1 in w 2 . All the results derived at n = 2 are easily extended. 
We leave the case where n is not integer an open question. 

As an example, we seek a pulse solitary wave in an unpre-stretched n = 3 
power-law dispersive solid. We find that the counterpart to (14. 2 p is 

3 27 4 

We look for the linearly polarized, pulse solitary wave corresponding to the 
case /3q 0, fli = 0, E = 0. To make a meaningful comparison with the 
n = 2 case, we perform the same changes of variable and function as in (14.51) . 
We find that the counterpart to (14.61) is 

<^-« 2 (r 2 +l^<H' (417) 

This differential equation can actually be solved using hyperbolic functions. 
Rather than present the details of that long resolution, we rapidly discuss the 
effect of having a "stiffer" power-law material — in the sense that n goes from 
2 to 3 while fi and b remain the same. When n = 2 the maximal amplitude of 
the wave is w(0) = 1, see (I4.6P ; when n = 3 the maximal amplitude is found 
by solving the bracketed biquadratic on the right hand-side of (I4.17p . An 
elementary comparison shows that the maximal amplitude at n = 3 is always 
larger than 1. Hence the "stiffening" of the material increases the amplitude 
of the pulse solitary wave. Clearly, the same conclusion can be reached for 
the compact wave corresponding to the case (3q = 0, (3\ ^ 0, E = 0. 
For the special case of fourth-order elasticity (I2.13p . we find that 

i 2 = -fi 2 /4, i 3 = 0, so that S(fi 2 ) = ^ 2 /4+z/ 4 fi 4 /16. (4.18) 
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Clearly, by comparing this equation with ( 14. ip . and by identifying z/ 4 with 
2b, we find the exact same results for the fourth-order elasticity theory of 
incompressible dispersive solids as we have for the n = 2 power-law solid. 
Hence in particular, transverse pulses and kinks with compact support are 
possible in a forth-order elastic solid with constitutive equations ( I2.13P and 

5 Concluding remark: 

A vector MKdV equation 

To conclude, we go back to the general governing equation (13. 9p . We do not 
specialize the constitutive relations for £ and a, but we perform a moving 
frame expansion with the new scales s = z — ct, r = et. 
We assume that W is of the form 

W = e 1/2 w, where w = 0(1). (5.1) 

Then Q — \W\ = e 1//2 |K;| and we expand the terms in (13.91) as 

(QW) ZZ = e^ 2 Q{^)w ss + e 3 / 2 Q'(0) {\w\ 2 w) ss + ..., 
pWtt = e 1/2 pc 2 w ss - 2e 3/2 pcw ST + 

{aW a ) zz = e l ' 2 c 2 a^))w ssss - 2e^ 2 ca{0)w SSST + ... (5.2) 

In order to recover the linear wave speed at the lowest order (here, e 1//2 ) given 
by Q(0) = pc 2 , we must assume that 

Q(0) = 0(1), and a(0) = 0(e) = ea (say), (5.3) 

where a is a constant of order 0(1). 
Then we find at the next order that 

0/(0) (\w\ 2 w) ss + c 2 a w ssss = -2pcw ST , (5.4) 

which we integrate once with respect to s to get the vectorial MKdV equation 
of Gorbacheva and Ostrosky [5], 

w T + q (\w\ 2 w) ss +pw sss = 0, (5.5) 

where here q = 0/(O)/(2pc) and p = ca /(2p). 

In this way we have clarified the status of the Gorbacheva and Ostro- 
vsky's beautiful results in the framework of the general theory of dispersive 
hyperelasticity. 
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